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Abstract. Given a digraph G = iVG, AG), an even factor M C AG is 
a subset of arcs tiiat decomposes into a collection of node-disjoint paths 
and even cycles. Even factors in digraphs were introduced by Geleen and 
Cunningham and generalize path matchings in undirected graphs. 
Finding an even factor of maximum cardinality in a general digraph is 
known to be NP-hard but for the class of odd-cycle symmetric digraphs 
the problem is polynomially solvable. So far, the only combinatorial al- 
gorithm known for this task is due to Pap; it has the running time of 
0(n*) (hereinafter n stands for the number of nodes in G). 
In this paper we present a novel sparse recovery technique and devise 
an 0(n^ log n)-time algorithm for finding a maximum cardinality even 
factor in an odd-cycle symmetric digraph. 



1 Introduction 

In [3] Cunningham and Geleen introduced the notion of independent path match- 
ings and investigated their connection to the separation algorithms for the 
matchable set polytope, which was previously studied by Balas and Pulley- 
blank [T]. Finding an independent path matching of maximum size was rec- 
ognized as an intriguing example of a graph-theoretic optimization problem that 
is difhcult to tackle by the purely combinatorial means. Two algorithms were 
given by Cunningham and Geleen: one relies on the ellipsoid method [4], and 
the other is based on deterministic evaluations of the Tutte matrix [5]. Later, 
a rather complicated combinatorial algorithm was proposed by Spille and Weis- 
mantel [13] . 

The notion of an even factor was introduced as a further generalization of 
path matchings in a manuscript of Cunningham and Geleen [5] (see also [3])- 
An even factor is a set of arcs that decomposes into a node-disjoint collection 
of simple paths and simple cycles of even lengths. Since cycles of length 2 are 
allowed, it is not difficult to see that finding a maximum matching in an undi- 
rected graph G reduces to computing a maximum even factor in the digraph 
obtained from G by replacing each edge with a pair of oppositely directed arcs. 
On the other hand, no reduction from even factors to non-bipartite matchings 
is known. 
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Finding a maximum cardinality even factor is known to be NP-hard in general 
digraphs [3] . For the class of weakly symmetric digraphs a min-max relation was 
established by Cunningham and Geleen |B] and by Pap and Szego [T2]. Later 
it was noted by Pap [10 that these arguments hold for a slightly broader class 
of odd- cycle symmetric digraphs. Takazawa and Kobayashi |16| pointed out a 
relation between even factors and jump systems and showed that the requirement 
for a digraph to be odd-cycle symmetric is natural, in a sense. 

The question of finding a combinatorial solution to the maximum even factor 
problem in an odd-cycle symmetric digraph stood open for quite a while until 
Pap gave a direct 0{n'^)-time algorithm |10lll| . His method can be slightly sped 
up to 0{n^{m + n\ogn)), as explained in Section |31 (Hereinafter n stands for 
the number of nodes in G and m denotes the number of arcs.) To compare: the 
classical algorithm of Micali and Vazirani for finding a maximum non-bipartite 
matching, which is a special case of the maximum even factor problem, runs in 
0{mn}^'^) time [5]. It is tempting to design a faster algorithm for the maximum 
even factor problem by applying the ideas developed for matchings (e.g. blocking 
augmentation [7]). There are, however, certain complications making even the 
bound of 0{mn) nontrivial. 

To explain the nature of these difficulties let us briefly review Pap's approach 
(a more detailed exposition will be given in Section [3]). It resembles Edmonds' 
non-bipartite matching algorithm and executes a series of iterations each trying 
to increase the size of the current even factor M by one. At each such itera- 
tion, a search for an augmenting path P is performed. If no such path is found 
then M is maximum. Otherwise, the algorithm tries to apply P to M. If no odd 
cycle appears after the augmentation then the iteration completes. Otherwise, a 
certain contracting reduction is applied to G and M . 

Hence, each iteration consists of phases and the number of nodes in the cur- 
rent digraph decreases with each phase. Totally there are 0{n) iterations and 
0{n) phases during each iteration, which is quite similar to the usual blossom- 
shrinking method. The difference is that during a phase the reduction may com- 
pletely change the alternating reachability structure so the next phase is forced 
to start looking for P from scratch. (Compare this with Edmonds' algorithm 
were a blossom contraction changes the alternating forest in a predicable and 
consistent way thus allowing this forest to be reused over the phases.) 

In this paper we present a novel 0{n^ log n)-time algorithm for solving the 
maximum even factor problem. It is based on Pap's method but grows the al- 
ternating forest in a more careful fashion. When a contraction is made in the 
current digraph the forest gets destroyed. However, we are able to restore it by 
running a sparse recovery procedure that carries out a reachability search in 
a specially crafted digraph with 0{n) arcs in O(nlogn) time (where the logn 
factor comes from manipulations with data structures). 

Our method seems applicable to a large variety of related problems. In par- 
ticular, the 0(mn'^)-time algorithm of Takazawa [2] solves the weighted even 
factor problem in 0{mn^) time and also involves recomputing the alternating 
forest from scratch on each phase. 



Similar effects of reachability failure are known to occur in the maximum 
C4-free 2- factor problem [TT]. Also, adding matroidal structures leads to the 
maximum independent even factor problem, which is solvable by the methods 
similar to the discussed above, see [S]. All these problems can benefit from the 
sparse recovery technique. Due to the lack of space we omit the details on these 
extensions. 



2 Preliminaries 



We employ some standard graph-theoretic notation throughout the paper. For 
an undirected graph G, we denote its sets of nodes and edges by VG and EG, 
respectively. For a directed graph, we speak of arcs rather than edges and denote 
the arc set of G by AG. A similar notation is used for paths, trees, and etc. We 
allow parallel edges and arcs but not loops. As long as this leads to no confusion, 
an arc from u to f is denoted by (u, v). 

A path or a cycle is called even (respectively odd) if is consists of an even 
(respectively odd) number of arcs or edges. For a digraph G a path-cycle matching 
is a subset of arcs M that is a union of node-disjoint simple paths and cycles in 
G. When M contains no odd cycle it is called an even factor. The size of M is 
its cardinality and the maximum even factor problem prompts for constructing 
an even factor of maximum size. 

An arc {u, v) in a digraph G is called symmetric if (v, u) is also present in G. 
Following the terminology from [10], we call G odd-cycle symmetric (respectively 
weakly symmetric) if for each odd (respectively any) cycle C all the arcs of C are 
symmetric. As already noted in Section [1] the maximum even factor problem is 
NP-hard in general but is tractable for odd-cycle symmetric digraphs. 

Maximizing the size of an even factor M in a digraph G is equivalent to 
minimizing its deficiency def (G, M) := \VG\ — \M\. The minimum deficiency of 
an even factor in G is called the deficiency of G and is denoted by def(G'). 

For a digraph G and U C VG, the set of arcs entering (respectively leaving) U 
is denoted by Sq{U) and Sq^^{U). Also, we write ^g{U) to denote the set of arcs 
with both endpoints in U and G\U] to denote the subgraph of G induced by U , 
i.e. G\U] — [U, jg{U)). When, the digraph is clear from the context it is omitted 
from notation. 

To contract a set U C VG in a digraph G means to replace nodes in J7 by a 
single complex node. The arcs in 7(V^G — U) are not affected, arcs in j{U) are 
dropped, and the arcs in S™{U) (respectively (5°"'(C/)) are redirected so as to 
enter (respectively leave) the complex node. The resulting graph is denoted by 
G/U. We identify the arcs in G/U with their pre-images in G. Note that G/U 
may contain multiple parallel arcs but not loops. If G' is obtained from G by an 
arbitrary series of contractions then G' = G/Ui/ . . . /Uk for a certain family of 
disjoint subsets Ui, . . . ,Uk C VG (called the maximum contracted sets). 



Algorithm 1 Simple- Augment(G, M) 



1: Search for an augmenting path P in G{M) 
2: if P docs not exist then 
3: return null 

4: else if P exists and is feasible then 

5: return M A A{P) 

6: else {P exists but is not feasible} 

7: Put M^ A A{P^) for i = 0, . . . , fc + 1 

8: Find an index i such that Mi is an even factor while Mj+i is not 

9: Find the unique odd cycle C in Mj+i 
10: G' ^ G/C, M' ^ m//C 
11: m' Simple- Augment(G', M') 
12: if M' = NULL then {M' is maximum in G'} 
13: return NULL 

14: else {M' is augmented in G' to a larger even factor M } 

15: Undo the contractions and transform M to an even factor M"*" in G 

16: return M+ 

17: end if 

18: end if 



3 Pap's Algorithm 



Consider an odd-cycle symmetric digraph G. The algorithm for finding a maxi- 
mum even factor in G follows the standard scheme of cardinality augmentation. 
Namely, we initially start with the empty even factor M and execute a series of 
iterations each aiming to increase \M\ by one. Iterations call Simple- Augment 
routine that, given an odd-cycle symmetric digraph G and an even factor M in 
G either returns a larger even factor or null indicating that the maximum 
size is reached. 



3.1 Augmentations 

Let us temporarily allow odd cycles and focus on path-cycle matchings in G. The 
latter are easily characterized as follows. Construct two disjoint copies of VG: 

:= {v^ \ ve VG] and := {v"^ \ v e VG}. For each arc a = {u,v) G AG 
add the edge {^1,^2} {corresponding to a). £)enote the resulting undirected 
bipartite graph by G. 

Clearly, a path-cycle matching M in G is characterized by the following 
properties: for each node v € VG, M has atjnost one arc entering v and also at 
most one arcjeaving v. Translating this to G one readily sees that M generates 
a matching M in G. Moreover, this correspondence between matchings in G and 
path-cycle matchings in G is one-to-one. A node (respectively m^) in G not 
covered by M is called a source (a sink, respectively). 



(a) 



(b) 



Fig. 1. Preparing for a contraction. Subfigure (a): the arcs of M are bold and 
the grayed arcs correspond to path Pi+i. Subfigure (b): path Pi is applied, the 
arcs of Mi are bold, and the grayed arcs indicate the remaining part of Pi+i- 

Given a digraph G and a path-cycle matching M in G we turn G into a 
digraph by directing the edges {y} , v^} corresponding to arcs (m, w) S M 

from to and the other edges from to . 

Definition 1 A simple path in starting in a source node is called alter- 

nating. An alternating path ending in a sink node is called augmenting. 

For an alternating path P let A{P) denote the set of arcs in G corresponding to 
the arcs of P in Hereinafter A l\ B denotes the symmetric difference of 

sets A and B. The next statements are well-known. 

Claim 1 If~d(M) admits no augmenting path then M is a path- cycle matching 
of maximum size. 

Claim 2 If P is an augmenting (respectively an even alternating) path in 
then M' :— M A A{P) is path-cycle matching obeying \M'\ — \M\ 1 
(respectively \M'\ — \M\). 

The augmentation procedure (see Algorithm [T]) constructs cf (M) and 
searches for an augmenting path P there. In case no such path exists, the 
current even factor M is maximum by Claim [T] (even in the broader class of 
path-cycle matchings), hence the algorithm terminates. Next, consider the case 
when ^(M) contains an augmenting path P. Claim [2] indicates how a larger 
path-cycle matching M' can be formed from Af , however M' may contain an 
odd cycle. The next definition focuses on this issue. 

Definition 2 Let P he an augmenting or an even alternating path. Then P is 
called feasible if M' := M A A[P) is again an even factor. 

If P is feasible then Simple- Augment exits with the updated even factor 
M A A{P). Consider the contrary, i.e. P is not feasible. Clearly, P is odd, say it 
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Fig. 2. Some cases appearing in Lemma [TJ Digraph G and even factor N (bold 
arcs) are depicted. 

consists of 2/c + 1 arcs. Construct a series of even alternating paths Pq, . . . ,Pk 
where Pi is formed by taking the first 2i arcs of P (0 < i < fc). Also, put 
Pfc+i := P and M, := M A A{Pi) (0 < i < fc + 1). 

Then there exists an index i [Q < i < k) such that Pi is feasible while 
Pi+i is not feasible. In other words, Mi is an even factor obeying def (G, Mi) = 
def (G,M) and M^+i contains an odd cycle. Since Mi and M^+i differ by at 
most two arcs, it can be easily showed that an odd cycle in Af^+i, call it C, 
is unique (see [10]). Moreover, Mi jits G, that is, \Mi<^AC\ = \VC\ - 1 and 
5°"*(yG) n Mi = 0. See Fig. [D for an example. 

It turns out that when an even factor fits an odd cycle then a certain 
optimality-preserving contraction is possible. As long as no confusion is pos- 
sible, for a digraph H and a cycle K we abbreviate H/VK to H/K. Also, for 
X C AH we write X/K to denote X \ "fniVK). 

Claim 3 (Pap jlOj ) Let K he an odd cycle in H and N be an even factor that 
fits K. Put H' := H/K and N' :~ N/K. Then H' is an odd- cycle symmetric 
digraph and N' is an even factor in H' . Moreover, if N' is maximum in H' then 
N is maximum in H . 

Note that M is maximum in G if and only if Mi is maximum in G. Simple- 
Augment contracts G in G. Let G' := G/G and M' := M^/G denote the 
resulting digraph and the even factor. To check if M' is maximum in G' a recur- 
sive call Simple- Augment(G', M') is made. If null is returned then M' is a 
maximum even factor in G', which implies by Claim [3] that the initial even fac- 
tor M was maximum in G. In this case Simple- Augment terminates returning 

NULL. 

Suppose that the recursive call was able to augment M' to a larger even 
factor M in G'. The algorithm transforms M to an even factor in G with the 
help of the following statement from [I^ (see Fig. [5] for examples): 

Lemma 1. Let K be an odd cycle in an odd- cycle symmetric digraph H . Put 
H' := H/K and let N' be an even factor in H' . Then there exists an even 
factor N in H obeying def(H,iV) = <:M{H',N'). 



3.2 Complexity 

There arc 0{ti) iterations each consisting of 0{n) phases. To bound the com- 
plexity of a single phase note that it takes 0{m) time to find an augmenting 
path P (if the latter exists). We may construct all the path-cycle matchings 
Mq, . . . , Mk+i and decompose each of them into node-disjoint paths and cycles 
in O(n^) time. Hence, finding the index i and the cycle C takes O(n^) time. Con- 
tracting C in G takes 0(m) time. (We have already spent 0{m) time looking 
for P, so it is feasible to spend another 0{m) time to construct the new digraph 
G' = G/C explicitly.) An obvious bookkeeping allows to undo all the contrac- 
tions performed during the iteration and transform the final even factor in the 
contracted digraph an even factor the initial digraph in 0{m) time. Totally, the 
algorithm runs in 0{'nf') time. 

The above bound can be slightly improved as follows. Note that the algorithm 
needs an arbitrary index i such that Mj is an even factor and Mj+i is not, i.e. 
i is not required to be minimum. Hence, we may carry out a binary search over 
the range [0, fc + 1]. At each step we keep a pair of indices (Z, r) such that M; is 
an even factor while Mr is not. Replacing the current segment [/, r] by a twice 
smaller one takes 0(n) time and requires constructing and testing a single path- 
cycle matching Mt where t := [(/ + r)/2\. This way, the 0{n^) term reduces 
to O(nlogn) and the total running time of becomes 0{n^{m + nlogn)). The 
ultimate goal of this paper is to get rid of the 0(m) term. 

4 A Faster Algorithm 
4.1 Augmentations 

The bottleneck of Simple- Augment are the augmenting path computations. 
To obtain an improvement we need better understanding of how these paths are 
calculated. Similarly to the usual path-finding algorithms we maintain a directed 
out-forest rooted at the source nodes. The nodes belonging to this forest are 
called T -reachable. At each step a new arc (u, v) leaving an J'-reachable node u 
is scanned and either gets added to J" (thus making v .F-reachable) or skipped 
because v is already J^-reachable. This process continues until a sink node is 
reached or no unscanned arcs remain in the digraph. 

Definition 3 Let G he a digraph and M he an even factor in G. An alternating 
forest J" for M is a directed out-forest in ^(M) such that: (i) the roots of T 

are all the source nodes in (ii) every path from a root of T to a leaf of 

J- is even. 

The intuition behind the suggested improvement is to grow carefully and avoid 
exploring infeasible alternating paths. 

Definition 4 An alternating forest T is called feasible if every even alternating 

or augmenting path in T is feasible. An alternating forest T is called complete 
if it contains no sink node and for each arc {u, v) in if u is T -reachable 

then so is v. 



We replace Simple- Augment by a more sophisticated recursive Fast- 
Augment procedure. It takes an odd-cycle symmetric digraph G, an even fac- 
tor M in G, and an additional flag named sparsify. The procedure returns a 
digraph G obtained from G by a number of contractions and an even factor M 
in G. Additionally, it may return an alternating forest J- for M in G. Exactly 
one of the following two cases applies: 

(1) def (G,M) = def(G, M) - 1 and T is undefined; 

(2) def (G, Af ) — def (G,M), M is maximum in G, M is maximum in G, and 

is a proper complete alternating forest for M in G. 

Assuming the above properties are true, let us explain how Fast- Augment 
can be used to perform a single augmenting iteration. Given a current even 
factor M in G the algorithm calls Fast-Augment(G, M, true) and examines 
the result. If def(G, M) — def (G, M) then by Af is a maximum even factor 
in G, the algorithm stops. (Note that the forest which is also returned by 
Fast- Augment, is not used here. This forest is needed due to the recursive 
nature of Fast- Augment.) Otherwise def(G,M) = def(G,M) - 1 by (U); this 
case will be referred to as a breakthrough. Applying Lemma [U M is transformed 
to an even factor M+ in G such that def(G, A/+) = def (G,M) = def(G, M) - 1. 
This completes the current iteration. 

Clearly, the algorithm constructs a maximum even factor correctly provided 
that Fast- Augment obeys the contract. Let us focus on the latter procedure. It 
starts growing a feasible alternating forest J" rooted at the source nodes (line [J). 
During the course of the execution, Fast- Augment scans the arcs of G in a 
certain order. For each node u in G we keep the list L(u) of all unscanned arcs 
leaving u. The following invariant is maintained: 

(3) if a = (m, v) is a scanned arc then either a G Af or both and are 
J^- reachable. 

Consider an J^-reachable node u^. To enumerate the arcs leaving in 
we fetch an unscanned arc a = {u,v) from L{u). If a G Af or is J-"-reachable 
then a is skipped and another arc is fetched. (In the former case a does not 
generate an arc leaving in C^(Af), in the latter case is already J^-reachable 
so a can be made scanned according to ([3]).) 

Otherwise, consider the arc oi :— (m^,w^) in Z^(Af) and let Pq denote the 
even alteranting feasible path from a root of 7^ to u^. Note that each node in 



(for X G VG) is either a sink or has the unique arc leaving it. A single step 



occurs when v"^ is a sink (lines fTOlTTSl) . The algorithm constructs an augmenting 
path Pi — Pq o ai leading to u^. (Here Li o L2 stands for the concatenation of 
Li and L2-) If Pi is feasible, the current even factor gets augmented according 
to Claim[2]and Fast- Augment terminates. Otherwise, forest growing stops and 
the algorithm proceeds to line [52] to deal with a contraction. 

A double step is executed when v'^ is not a sink flines fT5l- [20|) . To keep the 
leafs of J" on even distances from roots, T it must be extended by adding pairs 




Algorithm 2 Fast-Augment(G, M, sparsify) 



1: Initialize forest J- 

2: for all unscanned arcs o = (w, v) such that v} € VF do 

3: Mark a as scanned 

4: if a G M or i)^ G VT then 

5: continue for {to line 2} 

6: end if 

7: ai<^(w\v^) 

8: Let Po be the even alternating path to in {Po is feasible} 

9: if is a sink then {single step} 

10: P\ <= Pqo ai {Pi is augmenting} 

11: if Pi is feasible then 

12: return (G, M A yl(Pi), NULL) 

13: end if 

14: else {double step} 

15: Let a2 = {v^,w^) be the unique arc leaving {w^ ^ ^•^} 

16: Pi ^ Po o ai o a2 {Pi is even alternating} 

17: if Pi is feasible then 

18: Add nodes and and arcs ai, a2 to T 

19: continue for {to line 2} 

20: end if 

21: end if 

22: Mo^M AA{Po),Mi^M AA{Pi) 

23: Find a unique odd cycle C in Mi 

24: G' <S= G/C, M' <^ Mq/C 

25: if sparsify = false then 

26: return Fast- Augment(G', M', false) 

27: end if 

28: Construct the digraph H' 

29: (H',Ji',T) <;= FAST-AUGMENT(fl"', M' , FALSE) 

30: Compare VH and VG: let Zi, . . . , Zk be the maximum contracted sets and 

zi, . . . , Zk be the corresponding complex nodes in H 

31: G' ^G/Zi/_.../Zk 

32: if dcf(G', m') < def(G', M') then 

33: return (G',M', null) 

34: end if _ 

35: Unscan the arcs in G that belong to M and the arcs that enter zi, . . . ,Zk 

36: G ^G' , M ^M' , 7 

37: end for 

38: return {G,M,F) 



of arcs. Namely, there is a unique arc leaving in G(M), say a2 = (u^, w^) 
(evidently {w,v) G M). Moreover, is not a source node and is the 

only arc entering . Hence, is not J^-rcachablc. If Pi := Pooaioa2 is feasible 
then ai and a2 are added to thus making and J^-reachable. Otherwise, 
a contraction is necessary. 



Now we explain how the algorithm deals with contractions at line [221 One 
has an even alternating feasible path Pq and an infeasible augmenting or even 
alternating path Pi (obtained by extending Pq by one or two arcs). Put Mq := 
M A A{Po) and Mi := M A A{Pi). Let C denote the unique odd cycle in Mi. 
Put G" := G/C, M' := Mn/C. If sparsify = false then Fast- Augment acts 
similar to Simple- Augment, namely, it makes a recursive call passing G' and 
Af as an input and, hence, restarting the whole path computation. 

Next, suppose sparsify = true. In this case the algorithm tries to recover 
some proper alternating forest for the contracted digraph G' and the updated 
even factor M' . To accomplish this, a sparse digraph H' is constructed and 
Fast- Augment is recursively called for it (with sparsify = false). The latter 
nested call may obtain a breakthrough, that is, find an even factor of smaller 
deficiency. In this case, the outer call terminates immediately. Otherwise, the 
nested call returns a complete proper alternating forest for an even factor 
M in a digraph H (obtained from H' by contractions). This forest is used by 
the outer call to continue the path-searching process. It turns out that almost 
all of the arcs that were earlier fetched by the outer call need no additional 
processing and may remain scanned w.r.t. the new, recovered forest. This way, 
the algorithm amortizes arc scans during the outer call. 

More formally, the algorithm constructs H' as follows. First, take the node 
set of G, add all the arcs of M and all the arcs (u, v) G AG such that {u^,v^) is 
present in J". Denote the resulting digraph by H. We need to ensure that H is 
odd-cycle symmetric: if some arc {u, v) in already added to H and the reverse 
arc {v,u) exists in G then add {v,u) to H. Next, put H' :— H/C. Note that H' 
is a sparse spanning subgraph of G' (i.e. VH' = VG\ \AH'\ ~ 0{n)) and M' is 
an even factor in H' . 

The algorithm makes a recursive call Fast- Augment (i?', M', false). Let 
H and M be the resulting digraph and the even factor, respectively. Compare 
the node sets of H and H . Clearly, H is obtained from H hy a. number of 
contractions {G/C being one of them). Let Zi, . . . , Zfc be the maximum disjoint 
subsets of G such that h' = G/Zi/ . . . /Z^. Also, let zi, . . . , Zk be the composite 
nodes in H corresponding to Zi , . . . , . The algorithm applies these contrac- 
tions to G and constructs the digraph G := G / Zi/ . . . / Zk- Clearly M is an 
even factor in both g' and i?'. If def (i?',M') < def (i?', M') = def(G, M), then 
one has a breakthrough. Fast- Augment terminates yielding G and M . 

Otherwise, the recursive call in line[5S]also returns a proper complete forest T 
for H and M . Recall that some arcs in G are marked as scanned. Since we 
identify the arcs of G with their pre-images in G, one may speak of scanned 
arcs of G . The algorithm "unscans" certain arcs a = (u, v) G AG by adding 
them back to their corresponding lists L{u) to ensure ([3]). Namely, the arcs that 
belong to M and are present in G and the arcs that enter any of the complex 
nodes zi, . . . , in G are unscanned. After this, the algorithm puts G G , 
M := M , F :— F and proceeds with growing F (using the adjusted set of 
scanned arcs). 



Finally, if Fast- Augment has scanned all the arcs of G and is unable to 
reach a sink, the resulting forest J- is both proper and complete. In particular, 
by ([3]) no augmenting path for M exists. By Claim[3]this implies the maximality 
of M in G. The algorithm returns the current digraph G, the current (maximum) 
even factor M, and also the forest J^, which certifies the maximality of M. 

The correctness of Fast- Augment is evident except for the case when it 
tries to recover J' and alters the set of scanned arcs. One has to prove that ^ 
holds for the updated forest and the updated set of the scanned arcs. The proof 
of this statement is given in Section |^ 

4.2 Complexity 

We employ arc lists to represent digraphs. When a subset C/ in a digraph F is 
contracted we enumerate all arcs incident to U and update the lists accordingly. 
If a pair of parallel arcs appears after contraction, these arcs are merged, so all 
our digraphs remain simple. The above contraction of U takes 0(|F/^| • \U\) time. 
During FAST- Augment the sum of sizes of the contracted subsets telescopes to 
0{n), so graph contractions take 0{n?) time in total. The usual bookkeeping 
allows to undo contractions and recover a maximum even factor in the original 
digraph in 0{m) time. 

Consider an invocation Fast- AuGMENT(r, iV, false) and let us bound its 
complexity (including the recursive calls). The outer loop of the algorithm 
(lines EHTT]) enumerates the unscanned arcs. Since sparsify — false, each arc can 
be scanned at most once, so the bound of 0{\Ar\) for the number of arc scans 
follows. Using the appropriate data structures to represent even factors the reach- 
ability checks in lines [TT1 and [T71 can be carried out in 0(log It^-^l) time (see Sec- 
tion[B]for more details). Constructing Mo, Mi, and G takes 0{\Vr\) time. This 
way, Fast-Augment(/", N, false) takes 0((fc -I- 1) \ Ar\ log \Vr\) time, where k 
denotes the number of graph contractions performed during the invocation. 

Next, we focus on Fast- AuGMENT(_r, A^, true) call. Now one may need to 
perform more than \Ar\ arc scans since forest recovery may produce new un- 
scanned arcs (line [35]). Note that forest recovery totally occurs 0(|V^P|) times 
(since each such occurrence leads to a contraction). During each recovery M 
generates 0(|V^P|) unscanned arcs or, in total, 0{\Vr\^) such arcs for the du- 
ration of Fast- Augment. Also, each node Zi generates 0{\Vr\) unscanned 
arcs (recall that we merge parallel arcs and keep the current digraph sim- 
ple). The total number of these nodes processed during Fast- Augment is 
0{\Vr\) (since each such node corresponds to a contraction). Totally these 
nodes produce 0{\Vr\^) unscanned arcs. Hence, the total number of arc scans 

is oi\Ar\ + \vrf) = o{\vr\^). 

Each feasibility check costs 0(log \Vr\) time, or 0{\Vr\^ log \Vr\) in total. 
Finally, we must account for the time spent in the recursive invocations during 
FAST-AuGMENT(r, N, true). Each such invocation deals with a sparse digraph 
and hence takes 0{{k + 1) \Vr\ log \Vr\) time (where, as earlier, k denotes the 
number of contractions performed by the recursive invocation). Since the to- 



tal number of contractions is 0(|F-r|), the sum over all recursive invocations 
telescopes to 0{\Vrf log \Vr\). 

The total time bound for FAST- Augment (/^, N, true) (including the recur- 
sive calls) is also 0(|y/^|^ log |y/^|). Therefore a maximum even factor in an 
odd-cycle symmetric digraph can be found in 0(n^ log n) time, as claimed. 
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Appendix 



A Correctness 

In order to establish the correctness of the algorithm one needs to prove that 
once the forest T gets recovered at line [31] of Fast- Augment the set of scanned 
arc obeys property ([3]). The latter is equivalent to the following: 

Lemma 2. Consider line\^ in Algorithm\^ and let a — (u, v) be an arbitrary 
arc of G . Then either a is unscanned or a ^ M or both v} and are T- 
reachable. 

First, we shall need a more convenient characterization of reachable nodes. 
Consider a digraph F and a node v £ VF. Construct a new odd-cycle symmetric 
digraph F * v from F by adding a new node v' and an arc {v,v'). 

Lemma 3. For each F and v either dei{F * w) = def(/^) or def(I^ * = 
Aei{F) + l. 

Proof. 

Each even factor iV in is also an even factor in F * v, hence dcf(J^ * v) < 
def(r)-f 1. Also, for an even factor A^* in r*t; put N N*\{{v,v')}. Then, N 
is an even factor in 7^ obeying |A^| > |A^*| — 1. This implies def(-r) < dei{F*v), 
as required. □ 

Lemma 4. Let N be a maximum even factor in F and F be a feasible alter- 
nating forest for N . If is F -reachable then def(r' * w) = def(r'). 

Proof. 

Consider the even alternating path P from a root of T to w^. Since T is 
feasible, N' := N A A{P) is an even factor in F and no arc of N' leaves v' . Now 
N* :— N' U {(w, v')} is an even factor in F *v. Therefore, def (r * v) < def{F * 
V, N*) = dei{F, N) = dcf(r) implying def(r * w) = def(r) by Lemma[3l □ 

Lemma 5. Let N be a maximum even factor in F and F be a complete alter- 
nating forest for N. If is not T-reachable then def(-r * v) = dei{F) -\- 1 and 
N is a maximum even factor in F * v. 

Proof. 

The alternating forest T* obtained from T by adding a new root node {v'Y 
is complete for F * v and N . Hence, N is maximum in both F and F * v and 
de{(F*v) =def(r). □ 

Consider a node v 6 VG. If w € VG — [Zi U . . . U Z^) then we say that v survives 
the contractions. These nodes are both present in G and G . 



Lemma 6. Suppose that a node v e VG survives the contractions and is 
J- -reachable. Then is J- -reachable. Also, the nodes . . . , are T -reachable. 

Proof. 

The transformation oi H to H and of M to M may be viewed as follows: 

iH,M) {H,Mo) {H\M') (h\m') 
(4) II II II II 

(r",ivO) ^ (r",7vO) ^ {r\N^) {r\N^) ^ ... ^ {r',N') 

Here F'^, . . . are odd-cycle symmetric digraphs, N^, Nq are even factors in F"^ 
obeying |iV*| = |7V^|. Each fits some odd cycle C" in F' and T'+i ^ F'/C, 
iV^+i = N^/C\ 

Recall that the nested call to Fast- Augment in line [H] did not result into a 
breakthrough, so M is maximum in H and is a complete alternating forest. 

Let us prove the first claim of the lemma. Consider a node v surviving the 
contractions such that is J^-reachable in G. Suppose towards a contradiction 
that is not J^-reachable in H . By Lemma \Si M — A^" is a maximum even 
factor in 7? ^v^F^^v. Then, by Claim H N^'^ is a maximum even factor in 
F"'^ * V and hence so is N"'^ (by the equality of sizes). Proceeding this way 
in the backward direction we conclude that is a maximum even factor in 
F^ * V for alH = 0, . . . s. In particular, N'^ — M is maximum in F'~' = H . Since 
{v, v') ^ M, M is also maximum in H and def (i? *v,M)= def (F, M) + 1. This 
contradicts Lemma |4] and the fact that is J^- reachable. 

For the second claim, fix a node v = Zi and suppose that is not J-"- 
reachable in H . Consider the sequence of transformations (|4]) and suppose that 
V was created as a complex node in F^ while contracting an odd cycle C^~^ in 
F^~^ . As indicated above, is a maximum even factor in F^ * v. The latter, 
however, is false since A^"^ fits C^-i and, therefore, A^' has no arcs leaving v 
(cf. Fig. [IJb) for an example). Hence, A-' can be enlarged to A-' U {{v^ v')} — a 
contradiction. □ 

Lemma 7. Suppose that a node v £ VG survives the contractions and is 
T -reachable. Then is T-reachable. 

Proof. 

The node cannot be a source, hence it is reached by some arc (v} ^v'^^ 
in T , where is J^- reachable and a — (u, w) ^ M. Let ag — (uo, w) be the image 
of a under contractions. Note that a G AJF and ag G AH . By Lemma [6l is 
J^-reachable. Therefore, if oq ^ A/ then the completeness of T implies that v"^ 
is J^-reachable. It remains to consider the case ao G Af . The node Ug is not a 
source (since oq G M leaves uq) but is J^-reachable. In the auxiliary bipartite 
digraph u\ is entered by exactly one arc, namely (w^,mJ). Hence, must be 
J^-reachable, as required. □ 



Finally, we present the desired correctness proof. 



Proof of Lemma [2l 

Consider an arc a — {u, v) e AG . If a is not scanned then we are done. 
Otherwise let oq = (ug, vq) denote its pre-image in G. Here u — uq ii uq survives 
the contractions or u = Zi if uq £ Zi. Also, since all arcs entering the complex 
nodes zi , . . . , are unscanned (line I35p , vq survives the contractions and hence 
V = Vq. Since the algorithm only decreases the set of scanned arcs, ag must also 
be scanned in G. Clearly ao ^ M since all the arcs that belong to M and are 
present in G were as unscanned. Therefore, both uq and vq — v are J^-reachable 
in G by the invariant Applying Lemma [5] to uq and Lemma [7] to vq we see 
that both u and v must be J-'-reachable. The proof is now complete. □ 



B Feasibility Checks 

This section explains how path feasibility checks, which are performed by Fast- 
AuGMENT at lines [TT] and [ITl can be made efficient. That is, given an even 
factor N and a feasible even alternating path Pq we need to verify that an even 
alternating or an augmenting path Pi (obtained from Pq by appending one or 
two arcs) is feasible. We make use of a certain data structure V that maintains 
an even factor M A A{Po) as a collection of node-disjoint paths and cycles. The 
following operations are supported by V: 

— Insert(m,w): assuming that u is the end node of some path Pu in V and v 
is the start node of some path P„ in V, add the arc (u, v) thus linking P„ 
and Py or, in case P„ = Pv, turning this path into a cycle; 

— Remove(u, v): assuming that a — {u, v) is an arc belonging to some path or 
cycle in 2?, remove a thus splitting the path into two parts or turning the 
cycle into a path. 

— Is-Odd-Cycle(u, w): assuming that a = {u,v) is an arc belonging to some 
path or cycle in T), check if a belongs to an odd cycle. 

We make use of balanced search trees additionally augmented with Split and 
Concatenate operations (e.g. red-black trees, see |2|15| ) to represent paths and 
cycles in V. This way. Insert, Remove, and Is-Odd-Cycle take 0(log \Vr\) 
time each. Now checking if Pi for feasibility is done by calling Insert(u, v) and, 
in case of a double step, Remove(u>, v), and finally making Is-Odd-Cycle(u, v) 
request. If the latter indicates that Pi is not feasible, the changes in V are rolled 
back. 

During a Fast-Augment(P, N, false) call we grow J" in a depth-first fash- 
ion and maintain the structure 2? corresponding to the current J^-reachable 
node (i.e. T) keeps the decomposition of A A{P) where P is the path 
in J" from a root to u^). When T gets extended by arcs {u^,v'^) and {v'^,w^), 

becomes the new current node and V is updated accordingly by the above 
Insert(u,?;) and ILeuove{w,v) calls. When the algorithm backtracks from 
vj^ to u^, changes in V are reverted. This way, each feasibility check costs 
O (log I VP I) time 



Next, consider a FAST-AuGMENT(i", A^, true) call. The above time bound 
of 0(log |V^i^|) per check is only valid if we grow T from scratch. However, the 
algorithm also reuses the forest that is returned by the nested Fast- Augment 
call in line This incurs an overhead of 0(|yi^| log per forest recovery 

(this additional time is needed to traverse the arcs that are present in the recov- 
ered forest T and update T) accordingly). There are O(IT^r'l) forest recoveries 
during the call so the total overhead is 0(| log |y/^|) time. This does not 
affect the time bound of Fast- Augment (r. A, true). 



